Dynamics of Weakly Localized Waves 
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We develop a transport theory to describe the dynamics of (weakly) localized waves in a quasi-lD 
tube geometry both in reflection and in transmission. We compare our results to recent experiments 
with microwaves, and to other theories such as random matrix theory and supersymmetric theory. 



Localization of waves has always been among the most 
difficult yet most fascinating topics in the study of wave 
propagation in disordered media. The first studies dealt 
with infinite media, showing that localization is always 
achieved in ID, but that a minimum amount of disor- 
der is required in dimensions larger than 2 In 3D 
the critical point is estimated by the Ioffe-Regel crite- 
rion k£ rs 1, with k the wavenumber and £ the mean 
free path of the waves at a specified frequency @ ■ Later 
studies have considered localization in open media, 
and emphasized the 'leakage' through the boundaries — 
quantified by the conductance — as the basic localiza- 
tion parameter. The Thouless criterion 4] states that 
'leaky', extended states become localized when the 'di- 
mensionless conductance' g = G/(2e 2 /h) is of order one. 
Most recent studies, both in theory [S| and experiment 
0, have emphasized the giant fluctuations in transmis- 
sion coefficients in the regime g < 1, confirming the fun- 
damental importance of the Thouless conductance for all 
localization phenomena. In the diffuse regime (j > 1), 
apart from a factor of order unity, the dimensionless con- 
ductance g can be expressed as the ratio of the inverse 
microscopic level spacing, called the Heisenberg time tu, 
and the Thouless time t d = (L + 2z ) 2 /tt 2 Db (with Db 
the diffusion constant, L the size of the medium, and 
zq ~ £ accounting for internal reflection). 

A theory for 'all localization' does not exist. Important 
elements should be its capability to describe the tran- 
sition from the diffuse to the localized regime, notably 
with regard to leakage and dynamics, and in all dimen- 
sions, and its flexibility to experimental details, such as 
internal reflection, anisotropic scattering and absorption. 
A very complete localization theory is random matrix 
theory It can describe the transition from weak to 
strong localization, scaling, absorption, and fluctuations, 
and recently also dynamics but applies only for low- 
dimensional systems. Supersymmetric theory 0, [t| has 
also a very general range of applicability but does not 
always give the necessary physical insight to guide ex- 
periments. Finally, the self-consistent theory for local- 
ization 01 holds in all dimensions, is able to describe 
critical behavior around the mobility edge and has 
a clear generalization for dynamical problems. Its major 
disadvantages are that it applies only to the field correla- 
tion function and not to higher moment statistics, and its 



failure in the case of broken time-reversal invariance. It 
is valid on length scales larger than the mean free path, 
and — in the time domain — for times less than the 
Heisenberg time. 

Leakage effects can be studied from the 'leakage func- 
tion' (LF) Pt,r{ol) defined from the ensemble-averaged, 
time-dependent transmission (reflection) lT,R(t) accord- 
ing to, 



OO 

lT,R{t) = J dacxp(-at) P T ^ R (a) 



(1) 



Supersymmetric theories 0, Il2| have predicted strongly 
non-exponential decay in transmission, even for weakly 
localized waves (g ^S> 1) and in quasi- ID typically of the 
kind 'cxp[— g ln 2 (i/i# )]' beyond the Heisenberg time. It 
is in this regime that a modal picture is appropriate, like 
in chaotic cavities and that Pt,r{q) can be argued 
to equal the genuine distribution of resonant widths P(T) 
of the modes 0] at small T. P(T) has a log-normal be- 
havior at very small T attributed to 'prelocalized' modes 
0, which have become a central issue in the study of 
random lasers fl5|. 

For times smaller than the Heisenberg time tu super- 
symmetric theory predicts the transmission to decay like 
exp[— t/to + (I / gTT 2 )t 2 / t 2 D ] 0. This would imply a nar- 
row Gaussian distribution for Pr(a), centered around the 
average Thouless leakage 1/io with width ~ l/yfg~- A re- 
cent numerical simulation of wave dynamics in 2D disor- 
dered media has shown a similar, roughly quadratic 
increase of the logarithm of intensity. Chabanov, Zhang, 
and Genack recently studied weakly localized mi- 
crowaves in quasi-lD at times scales up to the Heisen- 
berg time, and observed a non-exponential transmission 
with time of the same type. Another interesting report 
— coming from random matrix theory |isj . and first re- 
ported for purely ID systems — is the l/t 2 reflec- 
tion coefficient for the semi-infinite quasi-lD tube, rather 
than the familiar l/t 3 / 2 decay expected from diffusion 
theory. This implies that in reflection Pr{o) cx a for 
small a. 

Transport theory ought to be valid for times less than 
the Heisenberg time, beyond which a modal picture takes 
over. The recent developments in theory and experiment 
call for a transport theory for the dynamics of (weakly) 
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localized waves, and notably for the leakage functions 
Pt,r(oz) denned in Eq. (JTJ. This is the subject of the 
present Letter. We will show that these functions are 
broadened by interference effects in a way compatible 
with observations and supersymmetric theory. We em- 
phasize that for a larger than the inverse Heisenberg time 
(the typical level spacing) the equivalence between the 
leakage function Pt,r(oc) and the resonant width distri- 
bution P(r) is not established, and that Pt,r{o) some- 
times takes negative values. 

Constructive interferences can be included into trans- 
port theory using the self-consistent theory of localiza- 
tion. In finite, open media this requires the appearance 
of a dynamical, spatially dependent diffusion constant 
D(r, 51) j2ij, which can explain the observed rounding of 
coherent backscattering of light near the mobility edge 
[20I l2l| , as well as the non-Ohmic transmission (23 • We 
will here study the dynamics. Given a short release of 
energy at the source at t — 0, the central observable is 
the flux of ensemble-averaged photon energy I(r, t) at 
position r and at time t, with Fourier transform I(r, 51) 
that we shall continue analytically in the whole complex 
plane. By causality is I(r, 51) an analytic function in the 
upper complex sheet Im 51 > 0. For positive times we can 
change the contour of the inverse Fourier transform with 
respect to frequency into the negative complex plane. If 
we assume that simple poles or branch cuts appear only 
along the negative imaginary axis, we find relation l[TJl 
with 
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FIG. 1: Time- dependent diffusion constant for wave trans- 
mission through a quasi- ID disordered waveguide. Theoreti- 
cal results (solid lines) are compared to experimental data of 
Ref. fl^l (dots). Satisfactory agreement between theory and 
experiment for times t below the Heisenberg time tn ~ goto 
is obtained by choosing go equal to 9 (sample A), 7.5 (B), 
and 4 (C) which is 14% to 33% larger than the experimen- 
tally reported values. Inset: The leakage function Pr(a) used 
to obtain the main plot. 



Pt,r{o.) = -ihm [I t ,r (H = -it 

cj.0 

- I t ,r (51 = -io. - e)] 



(2) 



In the normal diffuse regime only simple poles show up 
at 51 n = —in 2 ftp, and Pt,r(ck) equals an infinite sum of 
Dirac delta distributions. Purely localized modes would 
show up as a contribution 5(a) at zero leakage, but occur 
only in infinite or closed media. For an open quasi- ID 
system (N ^> 1 transverse modes, length L ^> £, classical 
diffusion constant Db — ve£/3, transport mean free path 
£ wavelength, and the energy transport velocity Ve) 
the basic equation is the ID dynamic diffusion equation 
for the intensity Green's function C(z, z', 51), 

[-m-d z D(z,n)d z ]c(z,z',n) = s(z-z'), (3) 

supplied by the self-consistency condition for the dy- 
namic diffusivity imposed by reciprocity [2oj |. 



(4) 



featuring the length scale £ = ^N£. At the boundaries 
z = 0, L we impose the usual radiative boundary condi- 
tions C =F z o [D(0/L,il)/D B )d z C = 0, where z ~ £ ac- 
counts for internal reflection. It.r is related to C through 
It,r{&) = TD(z = L/0, fL)d z C(z = L/0, z' = £, ft). 



The stationary problem (51 = 0) can be solved an- 
alytically by the substitution dr = dz/D(z, 0). This 
shows that for L ^S> £ the average transmission decays as 
exp(— L/£) which identifies £ as the localization length. 
The diffuse regime L <C £ has normal Ohmic transmis- 
sion with conductance g ~ go = t-N£/ (L + 2zo) ~ 2£/L. 
These results basicallyagree with the ones obtained from 
the DMPK equation [5j and supersymmetric theory [23j . 
Note that when L ^ £ it is important to discriminate be- 
tween go and the real conductance g which can be much 
smaller by localization effects. 

The solution for any complex- valued 51 has to be found 
numerically, by iteration. For go ^ 0.1, we found satisfy- 
ing and unique convergence for all 51 after 10-100 itera- 
tions. We have evaluated the leakage function by solving 
Eqs. (HUl for e = 10~ 9 /^ and e = l0~ w /t D and then 
using linear extrapolation to find the limit e J, 0. We have 
also carefully checked the absence of singularities away 
from the negative imaginary axis. The time-dependent 
transmission ir(£) was then obtained from Eq. Q). Ab- 
sorption can be added but this will just give rise to a 
trivial translation of the LF Pt(cc) to higher values for 
a. Following Chabanov, Zhang and Genack [J^ we shall 
interpret any non-exponential decay in terms of a time- 
dependent diffusion constant, in which case the transmis- 
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sion would decay as 



I T (t) ~ exp <^ - 



{L + 2z ) 2 



dt' D(t') 



(5) 



In Fig. Q] we have compared our calculation for D(t) to 
experimental results obtained for three different choices 
for the dimensionless conductance, corresponding to the 
samples A-C of Ref. g = 9 (A), g = 7.5 (B), 

and go = 4 (C). These values are slightly larger than can 
be estimated from the data of Ref. |ljj . Our transport 
theory describes the experimental results fairly well for 
all times below the Heisenberg time tn ~ goto- The 
inset of Fig. shows that the different branches of the 
leakage function Pt(o) achieve a finite width, though all 
with finite support. Note that the second branch has a 
negative value. We have fitted the first, positive branch 
to a Gaussian distribution with the same average and the 
same variance, and studied their variation with go. The 
Gaussian distribution leads to a linear decrease of D{t) 
shortly after the diffusion time tjj. Our findings can be 
summarized by the relation, 



D(t) A B t 

Db go go tD ' 



(6) 



with A = 0.15 and B = 0.20. Supersymmetric theory j^J 
gives B = 2/ ir 2 for orthogonal symmetry, in good agree- 
ment with our value, but makes no report of A. Yet, 
we have noticed that this term increases the agreement 
with experiment considerably. Our theory assigns no 
weight to Pt (a) for values smaller than a certain thresh- 
old a* w (l/i£>)(l — 0.8/,/on), in strong disagreement 
with supersymmetric theory |9j |l2j > which predicts a log- 
normal distribution for small a, caused by 'prelocalized' 
states that have localization lengths much smaller than 
the average localization length £. Our transport theory 
is not valid when a is small compared to the average level 
spacing. It is for this reason also that it can describe only 
the uninteresting short-time dynamics of transmission in 
the localized regime g < 1 where the Heisenberg time is 
smaller than the diffusion time. 

We will finally study the dynamics in reflection, and 
apply the same procedure to calculate the leakage func- 
tion Pr(o). For g > 1 wc find a series of clearly sep- 
arated branches, all positive in sharp contrast to trans- 
mission, and again with width ~ 1/ yjg. Their maxima 
typically vary as yfa which generates the typically diffuse 
1 / 1 3 / 2 tail in the time domain. The threshold leakage rate 
a* ~ 1/to causes an exponential decay at times beyond 
the diffusion time to- 

As g decreases the different branches of Pr(q.) start to 
join when go as 0.5. For go <C 1, the threshold leakage 
rate decreases exponentially with go: In a* ~ —1/go, and 
becomes rapidly very small, implying the disappearance 
of exponential decay. We will consider a waveguide of 
length L > f. We find that when L > 20£, P R (a) has 
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FIG. 2: Time-dependent reflection from a quasi- ID waveg- 
uide of length L>^: N = 20 (solid line) and N = 10 (dashed 
line). Dotted lines show the slopes and l/t 2 . Time has 

been normalized by the mean free time t 3 . The curves are 
obtained by Laplace-transforming Pr(o) shown in the inset 
(N = 20, L = 100^ — wavy solid line; N = 20, L -> oo - 
solid line; iV = 10, L —* oo — dashed line). 



converged to its asymptotic limit at L — > oo. In this 
limit, go = and tn = oo, and our theory applies at 
all times. The asymptotic Pr(o) roughly has a square- 
root behavior that is taken over by a linear slope for 
small values of a (see the inset of Fig. 0) . The linear law 
gives rise to the tail Iii(t) ~ l/t 2 in the time-domain, 
as can be seen in Fig. [3 This is consistent with the 
prediction of Titov and Beenakker using random matrix 
theory 18] . They have estimated the cross-over to occur 
at a time t ~ Af 2 i s (where t s is the mean free time), 
again consistent with our findings. We conclude that this 
interesting dynamical cross-over is well captured by the 
self-consistent transport theory, which, in contrast to the 
method of Ref. 01 > is n °t limited to the case of L ^> £ 
and can be applied to a waveguide of any length and at 
any time below the Heisenberg time tu- 
la conclusion, we have shown that the dynamics of 
(weak) localization both in transmission and in reflec- 
tion of a quasi- ID waveguide can be described by a self- 
consistent diffusion equation. This theory is not valid 
beyond the Heisenberg time, and other methods such as 
proposed by supersymmetric er-models have to be em- 
ployed. Stimulated by recent accurate time-resolved ex- 
periments of strongly disordered 3D materials close to the 
mobility edge |24|, a future challenge is the application 
of this theory to 2D and 3D systems. 
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